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Abstract 

We solve a hot twisted Eguchi-Kawai model with only timelike plaquettes in the 
deconfined phase, by computing the quadratic quantum fluctuations around the 
classical vacuum. The solution of the model has some novel features: the eigenvalues 
of the time-like link variable are separated in L bunches, if L is the number of 
links of the original lattice in the time direction, and each bunch obeys a Wigner 
semicircular distribution of eigenvalues. This solution becomes unstable at a critical 
value of the coupling constant, where it is argued that a condensation of classical 
solutions takes place. This can be inferred by comparison with the heat-kernel 
model in the hamiltonian limit, and the related Douglas-Kazakov phase transition 
in QCD2. As a byproduct of our solution, we can reproduce the dependence of the 
coupling constant from the parameter describing the asymmetry of the lattice, in 
agreement with previous results by Karsch. 



1 Introduction 

The twisted one-plaquette models, or Eguchi-Kawai models @] , have been exten- 
sively studied in the last fifteen years for their remarkable property of satisfying, 
in the large limit, the same loop equations as the corresponding lattice gauge 
theories. The reduction of the space-time degrees of freedom to internal degrees 
of freedom is the most startling property of these models. Their solution would 
amount to finding the master field for lattice gauge theories, and to solving them 
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in the large N limit. Unfortunately, in spite of the reduction to one plaquette, the 
saddle point structure of the Eguchi-Kawai model in quite complicated, due to the 
effect of the twists, and a complete solution appears to be beyond the reach of our 
present techniques. 

In this paper we study a twisted one-plaquette model, which is much simpler 
than the complete twisted Eguchi-Kawai model, but whose solution can be explicitly 
found, at least in the weak coupling (or high temperature) regime, in terms of 
eigenvalue distribution. The model under consideration is defined by the partition 
function 



and in the large N limit it is equivalent, as shown later in the paper, to a c? + 1 
dimensional lattice gauge theory at finite temperature, with L links in the compact- 
ified time dimension, and with only the couplings induced by time-like plaquettes. 
Hence the action described by eq. (|l|) arises naturally in finite temperature lattice 
gauge theories as a zero order approximation in which space-like plaquettes are 
completely neglected. This approximation has already been discussed in previous 
papers Q and it will be reviewed in the following section. 

The main result of our paper consists in solving the model in the broken phase, 
which occurs for values of j3 larger than a critical value Per and is characterised by 
a non vanishing vacuum expectation value of the Polyakov loop V = TrV^^. The 
solution has some novel, interesting features, namely an eigenvalue distribution for 
the matrix V which consists of L separated bunches centered around the L^^ roots 
of unity. The eigenvalue distribution within each bunch is given, with very good 
approximation, by a semicircular Wigner distribution of calculable radius. The 
splitting of the eigenvalue distribution into several bunches is a new feature that 
makes the model in consideration interesting in itself and susceptible of applications 
in other branches of physics, possibly for instance condensed matter physics. 

The paper is organised as follows. In section ^ the twisted Eguchi-Kawai reduc- 
tion in 3 + 1 dimensions is considered in the framework of finite-temperature lattice 
gauge theories, and a scheme of dimensional reduction is proposed that allows to 
extend the analysis to c? -|- 1 dimensions. In section ^ we derive the effective model 
for the Polyakov loops in the deconfined phase, up to quadratic order in the fluc- 
tuations. In section ^ the solution of the effective model and the phase transition 
occurring at a critical value of the coupling are discussed for a model with arbitrary 
number L of lattice links in the time direction. In this context we also derive the 
rescaling of the coupling j3 with the parameter p, describing the asymmetry of the 
original lattice, and we find results in excellent agreement with the ones by Karsch 
in 0. In section ^ we present our conclusions. 



Z= DV l[DU^d 
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2 Dimensional reduction in finite temperature large 
N LGT 



2.1 Twisted reduction in 3 + 1 dimensions 

The symmetric twisted Eguchi-Kawai model (TEK ) in 3 + 1 dimensions is a one 
plaquette model which is equivalent, in the large limit, to a lattice SU(N) pure 
gauge theory defined in a periodic box of size L , with = L^. Its action is given 
byQ: 

5sTEK = iV/5 5]ReTr(e'-f/^t/,t/tf/t) . (2) 

In the large limit the size of the box goes to infinity in all directions. If we want 
to describe a finite temperature gauge theory we need a lattice whose size in the 
compactified time direction remains finite also in the large limit. This presents a 
problem in the TEK approach where the size of the box is linked to A^. The problem 
can be overcome in two different ways. We can define the theory in an asymmetric 
box, whose spatial extent goes to infinity in the large A^ limit while its extent in 
the time direction remains finite. Several of such boxes with the corresponding 
twists exist in the literature p, 0. On the other hand, we can keep a symmetric 
box of size L but with different lattice spacings a and a,- respectively in space and 
time directions. If T is the temperature, then 1/T is the size of compactified time 
dimension, given by: 

- = Lar . (3) 

Let us introduce the asymmetry parameter p = a/ar- Then from (j^) we have 

p = aTL . (4) 

From eq. (|^) it is clear that in order to keep T finite the parameter p has to go 
to infinity like L = N^^'^ in the large A^ limit. The introduction of different lattice 
spacings is equivalent to the introduction of different couplings [3s and [3t for space- 
like and time-like plaquettes respectively. The relation between these new couplings 
and the asymmetry parameter p is given by 

A(p)=p(Ayn.m + 4c.(p)) ; = l^Z!!l!!L±i^ , (5) 

where /3symm is the coupling corresponding to a symmetric lattice, namely to p = 1. 
The terms proportional to /5symm in (i) stem from the requirement that a Yang- 
Mills theory is reproduced in the naive or classical continuum limit. The functions 
Cr and Co- on the other hand represent the result of quantum corrections and have 



-'^In this formula can assume the values 0, 1,2,3. Afterwards we will always utilise these 
indices to enumerate the space dimensions only. 
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been evaluated by F. Karsch in This point will be discussed more in detail in 
the section |4^ . 

The resulting one-plaquette model in this approach is then described by the 
following partition function: 

^hTEK = J DV J f[ DU^ exp (S'hTEK) 

3 3 

^hTEK = N(3tip)J2^eTTie''^VU^V^Ul) + Nl3,ip) ^ ReTr (e'^f/^f/,?7tf/t) 

(6) 

where we have denoted with V the link variable in the time direction, and with 
Ufj, the ones in space directions. It is important to remember that only traces 
corresponding to closed loops in the original theory^ have physical meaning, like for 
instance the Polyakov loop which is given by 

p = i_Tr\/^ . (7) 

We shall not consider here in detail the models based on the use of asymmetric 
boxes, which have the advantage of having Ps = Pt = P but the disadvantage of far 
more complicated twists, at least for the space-like plaquettes. It should be noticed 
however that, in the approximation where space-like plaquettes are neglected, the 
model based on asymmetric boxes leads to the same action given in (|I]). 

2.2 The perturbative expansion in /Sg 

If we were able to perform the exact integration over the space-like link variables 
in eq. (||), we would obtain the exact effective action for V, and thus for 
the Polyakov loop, in the large limit. However this is too difficult, and we 
shall follow a perturbative approach consisting in an expansion of the partition 
function (j^) in powers of Ps- In the present paper we shall consider only the zeroth 
order of this expansion, namely the partition function (|l]) which corresponds to 
neglecting the space-like plaquettes altogether. An explicit calculation of the first 
non trivial order in P^ seems to be within reach of the techniques developed in the 
following section, but it will be left to future investigation. The extent to which 
the zero order approximation is reliable was discussed in ref. [0], where the same 
approximation was used to study, in the large limit, the effective action for the 
Polyakov loop in a lattice gauge theory with both a Wilson and a heat kernel action. 
The difference between the approach of the present paper and the one of ref. 0] 

^That is, the d + 1-dimensional LGT (defined on a box of size L) from which the model (||) [or 
the model (^] is obtained by eliminating the space-time dependence of the gaug e fie lds via the 



"twisted Eguchi-Kawai" reduction (see 0, and similar arguments in next section 2.3 ) 
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is that the reduction to a one plaquette model was obtained there, for the Wilson 
action, by an approximate Migdal-Kadanoff bond-moving scheme. Instead our 
starting point eq. is exact. 

The perturbative expansion in (3s was studied in even more details in ref. 0], 
where the expansion was carried out, for the SU{2) lattice gauge theory, up to the 
order (3l. The feature that emerges from ref.s 0, ^ is that the zeroth order approx- 
imation, although it correctly predicts the existence of a deconfinement transition 
for large temperature, fails to reproduce the deconfinement temperature in 3 -|- 1 
dimensions for lattice sizes0 nt in the time direction greater than 2. This is related 
to the fact that in the zeroth order approximation the critical value of j3t rescales 
linearly with rit rather than with the logarithmic law required by asymptotic free- 
dom. A correct prediction of the latter behaviour would presumably require an 
exact derivation, namely to all orders in (3s-, of the effective action for the Polyakov 
loop; however in the SU{2) theory the first non trivial order in (3s already provides 
a striking improvement in the scaling properties We expect the same to hap- 
pen in the large N case, but in this paper we shall restrict ourselves, as already 
mentioned, to the zeroth order approximation, the emphasis being more on the the 
matrix model (|l]) itself and its interesting features. It should be noticed also that, 
unlike the 3 + 1 dimensional case, in 2 + 1 dimensions the full theory has the same 
linear rescaling with nt as the (3s = approximation, and that there are indications 
that such approximation is in that case reliable also for large rit. 

2.3 Twisted reduction in d + 1 dimensions. Weak coupling 
vacuum 

We have so far considered the reduction to a one plaquette model of finite tempera- 
ture lattice gauge theory only in 3 + 1 dimensions. We want now to show here that 
the partition function given in (|I]) can be obtained as the one plaquette reduction 
of the large N finite temperature lattice gauge theory in d + 1 dimensions, in the 
zeroth order approximation (3s = 0. For even space-time dimensions we could prove 
this by using the existing twisted Eguchi-Kawai models, essentially in the same way 
as in 3 + 1 dimensions. However in odd space-time dimensions consistent twists have 
not been constructed, and we have therefore to rely on a different type of approach. 

We consider a lattice gauge theory defined on a lattice which is infinite in the 
space directions and of length L with periodic boundary conditions in the time 
direction. The strategy is to reduce the height of the lattice in the time direction 
to one by using the same type of approach as in the Eguchi-Kawai model, and then 
to eliminate the dependence from space coordinates by assuming that a constant 

■^The lattice size nt introduced here corresponds to the number of hnks in the time direction 
in a lattice with the same lattice spacing in the space and time directions. So it is related to the 
lattice size L of the Eguchi-Kawai reduction of eq.(^ by the relation nt = ^. 
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master field dominates in the large limit. Let us start with the action 



S{L) = NY^Y^^eTi [PtT.Vix,t)U^{x,t + l)V^{x + fi,t)Ul{x,t) 

X t=l \ /^=1 

d 

+ /5s E U^{x,t)U^{x + fL,t)Ul{x + i),t)Ul{x,t) 

Notice that in (||) we kept two independent couplings Ps and Pf Unlike the case of 
eq. d^), however, this is not necessary in order to be able to take the limit —>■ oo, 
because L and are now unrelated, and we will be free to choose (3s = Pt at our 
convenience. The naive prescription for the reduction of the time degrees of freedom 
would be 

V{x,t) ^V{x) , 

U^ix,t) ^U^ix) . (9) 

It is easy to show by standard methods (see for instance ) that the resulting 
action, which can be identified with 5'(1), leads to the same set of loop equations 
in the large N limit as the full S{L) theory, provided all loops which are closed in 
the reduced lattice (L = 1) but correspond to open loops in the original lattice have 
vanishing expectation value. Just as in the untwisted Eguchi-Kawai model, this 
would be granted by the symmetry 

V{x) ^e'^V{x) (10) 

which is satisfied by >S'(1). The trace along "open" lines is not invariant under the 
symmetry ([l0| ) as they do not contain the same number of V{x) and l^^(x) fields. So 



these contributions vanish unless the symmetry is broken. The symmetry however 
is actually broken in the weak coupling regime, that is also in the continuum limit, 
where V{x) is close to 1 (more generically to an element of Z^y) and the traces of 
open lines do not vanish. Consequently the reduction prescription must be endowed 
with twists, namely it must be of the type 

V{x,t) D{x,t)V{x)D\x,t) , 

U^ix,t) ^ Dix,t)U^ix)D\x,t) , (11) 
where D[x,t) is given by 

D{x,t) = iT)T.-^iToy , (12) 
with r and Tq traceless SU{N) matrices satisfying the 't Hooft algebra 

TTo = e'^"'ToT . (13) 
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In the last equation m is an integer number to be determined. By performing in 
® the replacement (|TTD and redefining the variables according to the substitution 
U^{x) — > Ufj,{x)T and V{x) — > V{x)To we obtain 



= Jl[[DV{x)l[DU,{x)]exp{SR) 



d 

n 

Sr = AiV^^Re e^'^Tr[?7^(f)\/(x + /i)t/t(f)\/t(j)] + 



d 



[3,NY.Y.^^'^AU,mU,{x + fi)Ul{x + P)Ul{x)] . (14) 

Notice that, unlike the twisted Eguchi-Kawai model, the twists are present in ([T^) 
only in the contributions from time-like plaquettes, as the reduction has been done 
only in the time direction. 

Consider now the loop equations for the reduced theory ([l^ . We already re- 
marked that as long as the symmetry (p!OD is unbroken the loop equations of the 
reduced theory coincide with the ones of (H). We will show now that in the twisted 
theory this is the case also in the weak coupling regime. Indeed in the extremely 
weak coupling the fields tend to their vacuum configurations, which in our twisted 
reduced theory are 

Ui{x) PN®lm , 

m 

v{x) ^ g^®i™ , (15) 

where 1^ is the m x m unit matrix and Pn_ and Qn_ are the usual building blocks 

m m 

for the twist eating configurations: 

{P)ab = 6a+i,b ; {Q)ab = Sabe'"-^'' , a,b = l,...N/m (16) 

with periodicity in the the indices a and b, namely a = — + 1 means a = 1. It is 
clear from the context that we have to restrict the values of N so that N/m is an 
integer. With these vacuum configurations the trace of open lines is proportional 
to Tr((5jv)*, where t is the difference between the number of V^'s and y^'s in the 

771 

trace, namely the difference between the time coordinate of the initial and the final 
point of the path in original unreduced lattice. It is elementary to see from eq. (|T6|) 
that 

Tr(giv)* = unless t = k— (17) 

m m 

for integer k. In the unreduced lattice the closed loops correspond to t = kL, so 
the comparison of the two equation determines m : 

N , , 

m = — . (18) 

LJ 
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Fig. 1: The weak coupling solution in the large-jj consists of L distributions, 
centered at the multiples of^. The width of the distributions grows as (3 decreases, 
until finally an instability is attained at a certain critical value [3^ ■ 



With the above replacement eq. (M) is, in the large N limit, an exact dimensional 
reduction of (j^) on an infinite (i- dimensional lattice. In order to obtain the matrix 
model (|I]) that we intend to study, two more steps are required. First, we assume 
that an af-independent master field dominates the functional integral so that the 
field variables in (0) can be replaced by constants; second, we neglect, according 
to the previous discussion, the contributions of the space-like plaquettes by setting 
j3s = 0; correspondingly, from now on we write simply (3 for jSt- For d = 3 the final 
result coincides with the one obtained directly from the hot twisted Eguchi-Kawai 
model, and the same can be shown to happen for any even space-time dimension. 
The derivation in the present section extends the validity of the reduction to odd 
space-time dimensions. 

The vacuum configuration given in ([1^) and (|T6|) consists for V{x) of L bunches 
of ^ degenerate eigenvalues at the L^^ roots of 1. This means that the open Polyakov 
loop coincides in this extreme weak coupling limit with the unit matrix'*. 

In the next section the quadratic quantum fiuctuations around this vacuum are 
considered and their contribution is calculated, resulting into a broadening of the 
eigenvalue distribution within each bunch and eventually, as Pt decreases, into a 
phase transition towards a uniform eigenvalue distribution. 

3 Derivation of the effective model 

Let us consider the partition function (|l]) of the twisted one-plaquette model. The 
integrations over the unitary matrices f/^ are all independent and the partition 
function can then be rewritten as 



DV 



DU e 



/3Af Re e' T Tr (y 1/ yt f/t ) 

''More generically, we could in (|l5|) replace Im = with an element of Zn . Correspondingly, 
V{x)^ would be an element of Zat, which is indeed the generic situation in the broken phase. 
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1=1 



(19) 



where in the last hne we have gauge-rotated the matrix U to diagonahze V: 

V ->diag(e'^o,...,e'^^-0 . (20) 

J^{0) is the Haar measure of SU(N) expressed in terms of integration over the 
eigenvalues e'^* and it is given by 

ne)=Y[Asm\^-^) . (21) 



We want now to solve the model (|T9D in the weak coupling limit, by expanding 
around the vacuum configuration in whichQ 

V = Ql®-^nil ; P = lTr\/^ = l . (22) 

In this vacuum the eigenvalues of V are organised in L bunches, each composed of 
^ identical values e^, a = 1, . . . L. Accordingly we write the eigenvalues of the 
matrix V in the weak coupling phase as 

27rci 

= Oa,a = -J^ + (l)a + ^a,a , (23) 

with a = l,...L,a = l,...^ and with Lpa^a constrained by Y.a H^a^a = 0. We have 
separated here the degrees of freedom 0a corresponding to the fluctuations of the 
centres of the bunches from the fluctuations (fa,a within each bunch. 

We will show that in the large N/ L limit the model is solved, for f3 bigger than 
a calculable critical value /?cr, by an eigenvalue distribution for V consisting of L 
separated distributions, all centered at a multiple of In the large L limit we 
will show that the L bunches of eigenvalues all have the same Wigner semicircular 
distribution. 



3.1 Integration over the space-like links 

The first step in order to solve the model is to perform explicitly the integration 
over the unitary matrix U in the expression (|19|) of the partition function. This will 
lead us to an "effective" action for the eigenvalues of the matrix V, which will be 
studied by means of matrix models techniques. 

^Here again we assume ^ to be an integer. This assumption is justified, as in the Twisted 
Eguchi-Kawai model, by the behef that the large N and the continuum limit (large L) are smooth. 
It would nevertheless be of some interest to study the case of iV = mL + r with m and r integers 
and < r < L. Although sub-leading in the large N expansion, the dynamics of the r eigenvalues 
that cannot accommodate into L identical bunches might reveal some interesting features, and it 
is likely to depend on the existence of common divisors of L and r. 
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Let us consider the slightly more general integral 

I{e,^)= [dU e^""^^^""'^''-^'^^^''^' , (24) 



which reduces to the one of ( ]T9| ) for ipj = dj — 

In spite of its analogy with the Itzykson-Zuber integral, the integration in (p^) 
cannot be done exactly (except for iV = 2) as the Duistermann-Heckmann theorem 
does not apply in this case. However in the large N limit it is possible to apply to 
it the saddle-point method. In this way Kogan et. al ||^, by taking into account, 
like in the Itzykson-Zuber integral, all extrema - both maxima and minima - of the 
exponent, obtained the following asymptotic formula: 

f DUe^'. = exp{£.A.m} ,,^o(h) . (25) 

P being a permutation of the indices {i}. By replacing in (EH) Aij with f3N cos{6i — 
ipj) and by using the trigonometrical identity 

cos(6'i - ipp{i)) + cos{6j - tpp{j)) - cos(6'i - ipp{j)) - cos{6j - ipp{i)) 

= 4 sin(^) cos( ''+'^"^7'"^"'^> ) (26) 

we would obtain the following expression for the large limit of I{6,iIj): 

N{N — 1) 

where C = (PN) 2 is an irrelevant overall factor, which will be discarded in 
what follows. 



In applying the saddle point method to the integral at the l.h.s. of eq. (25) one 
should be aware that the gaussian integration around the extrema would diverge 
as a maximum turns into a minimum of the integrand. These divergences are 
signalled by the zeros at the denominator of the r.h.s of eq. (0). The zeros of 
J^{0) and J'l'ip) however are cancelled by the zeros at the numerator when the sum 
over permutations is performed, in analogy to what happens in the Itzykson-Zuber 
integral. Instead the zeros true singularities of the r.h.s 

of (P?!). These singularities are not present in the original integral and they denote 



that the equation (pT]) cannot be applied in such circumstances unless one restricts 



the sum to a suitable subset of permutations. 



Let us go back to our original integral of eq. ([T9|), which is obtained from ([27|) 



with the substitution ipi —>■ 9i — The vacuum configuration of such a model was 
discussed in the previous section and it is given by eq. (|15|), with N/m = L. 

To consider just quantum fluctuations around this vacuum is equivalent to re- 
strict the sum over permutations at the r.h.s. of ([27| ) to the permutations that map 
each bunch of eigenvalues into the next bunch, namely: 

P '■ da,a 0a+l,Pa(a) ■ (28) 



10 



where Pa(«) corresponds to a permutation of the ^ eigenvalues of the bunch at 
The other permutations are exponentially depressed compared to the ones above; 
in fact for each eigenvalue mapped from the bunch a to the bunch a + s + 1 the 
exponential in (^Tf ) is depressed by a factor 

exp /?Ar(l - cos — ^) ~ exp -(3N—-^ . (29) 



L2 

In the next section we shall precisely examine the quantum fluctuations around the 
classical vacuum given in (13), and determine for what value of /3 such fluctuations 



make the vacuum unstable. In other words we shall neglect all extrema of the action 
other than the ones given in eq. (pSf), and also keep the center of the bunches fixed 
at their classical value, namely we shall fix = through all the calculation. 
Of course the question arises whether a condensation of the other extrema, or an 
instability in the dynamics of the bunches as a whole may cause a transition before 
the vacuum becomes unstable due to the quantum fluctuations. In the next section 
we shall give strong arguments, based on the comparison of the present model with 
QCD in two dimension on a cylinder, that at least for large L the value of (3 for 
which the vacuum becomes unstable in very close, if not coincident, with the one 
where the transition takes place. Finally, let us remark that when the sum over 
permutations in eq. ( P7| ) is restricted to the ones of (^) the result of the saddle 
point method is singularity free, as the argument of the cosine at the denominator 
is for those permutation close to zero. 

3.2 The effective quadratic model 

With the replacement ipi 6i — ^ the r.h.s. of eq. (^Tf ) becomes 

V ^ ^ ^ n„ cos(^i^^^^i4^^) 



(30) 



where the sum is now restricted to the permutations given in eq. (Pq). We replace 
9i with its classical value plus fluctuations, according to eq. (^), and expand 
each term up to quadratic order in the fluctuations. It is easy to see that upon 
this expansion the fluctuations 0^ (the "centres") and ipa,a decouple completely. 
As already discussed, we consider the solution where all the (pa vanish. In the 
remaining of this section we keep therefore track only of the terms containing (pa,a 
and we derive a quadratic effective model for them. 

Let us consider first the factor J'^{ff) appearing both in the integration volume 
and in (|30|) . Under the substitution (|23|) it becomes 



jr2(^) = n4sin2^^^^ n 4sin2 

ij ^ a,a<b,l3 



7r(a - b) ^ ipa^a - ipb,p ' 



11 



X 



^^P] sin 



7r{a -b) ^ (Pa,a - (Pb,f3 



sm 



L 

'Pa,a—'Pb,l3 



'Pa,a—'Pb,l3 



(31) 



where A'^{{pa) = Y[a<p{¥^a,a — '^a,i3)^ IS the usual Cauchy-Vandermonde determinant. 
The quadratic expansion of eq. (PT|) is 



1 AT L—l -| 1 

[HA^M X exp{-if (E-k + ^)E<. + 0(^^)} . (32) 

r, ^ ±J / — ]^ Sin 1^ /I ^, 



The finite sum appearing above can be exactly calculated (see App. ^ and it is 
given by 



L-l 



E 



1 sin^ f 



(33) 



The exponential and the denominator in ([30| ) become respectively 

27r 

^ COS(f^fc - &p(^k) + 

k 



exp{^iV/?^C0s(6'fc - 9p^k) + -p)}} ^ exp{/3A^^ COs(v9a,a - Va+l,Paia))} 
i k a.,a 

= exp{-pN{J2vl^ -Y.Va,aV-+^,PM) + 0{v% (34) 



and 



COS ■ 



exp|— ^ log cos [ 

a,a<b,l3 

■N 



(35) 



Putting all these results together, we can now write down the explicit expression of 
the action (pj]) expanded up to quadratic terms in the fiuctuations ip: 



(1-d) 



d-l 



exp 



12 



X 



f^(-l)'^(^'') exp{^(/5L-^);^^,,«^,+i,p„(„)})' . 

\{Pa} ^ ^ '^.^ / 



(36) 



The solution of this model is the subject of the next section. 
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4 Solution of the model 



4.1 Zt, invariant solution 



In solving the model (|36|) for an arbitrary value of L, we assume that the master field 
is translational invariant, namely invariant under the Z^, symmetry of the vacuum 
(p^). With this assumption all bunches of eigenvalues have the same distribution 
and we can set: 

Va,a ^ , (37) 

where is now of order 1 in the large L limit (recall that the fluctuations (pa^a 
are of order j^. The invariance under translations is at this stage an assumption, 
although a very reasonable one; it will be checked a posteriori in section ^]3| for an 
arbitrary value of L. 

With the position (|37|) the partition function (^) reduces to the product of the 
partition functions of L identical models: 



KM 



/'I' T _ J W 1 

nd^4A2(v.)] ~ exp{-n[/3Hrf-^]E(^.r} 
a=l ^'^ a 

, (^{-ly^P^ exp{n/?^, ^c^P.i^) }) 



(3J 



where we have defined 



H 



(39) 



and n = N/ L. The partition function of each individual model in the above equation 
coincides with the one of a (i- dimensional Kazakov-Migdal model [^] with quadratic 
potential. In fact from the Kazakov-Migdal partition function 



/ UUDU.ix) I llD^x) exp{-n/3^Tr5:[^$2(f) 



^<l>(f)f/^(f)$(x + /i)f/t(f)]} 

11=1 



(40) 



by integrating over the unitary matrices U^{x) one obtains in the mean-field ap- 
proximation ($(a;) independent of x) the expression Zkm in (P^), where {v^a} are 
the eigenvalues of the hermitean matrix $, and the mass m is: 



ra 



2d- 



d-1 



(41) 



The model (|40|) was solved in the large-n limit by Gross |TT|, who showed that the 
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eigenvalue distribution p{ip) is semicircular: 



2 

p{ip) = — ^-y^^ — for W\ < ■ 



(42) 

p{(f) = for r < \ip\ < IT 



with radius given by 



2 4(2d-l) , 
r = 7 ^ — r • (43) 



(5H\m'^{d - 1) + rf^m^ - 4(2d - 1 

The solution becomes unstable at the value of (3h for which the radius become 
complex, namely 

= 12(d-VM-l) • <^^' 



The corresponding value of (3 (recall the definition (|39[)) represents a lower bound 



for the true critical coupling /?cr at which a phase transition takes place. 

It should be noticed however that a different upper limit may be found by looking 
at the value of j3H for which the radius of the eigenvalue distribution of becomes 
vr, corresponding to a radius ^ for each bunch of eigenvalues. A phase transition 
occurs at this point as the bunches merge together and the semicircular distribution 
is not anymore a solution. 

A better understanding of the phase transition can be obtained by noticing that 
for large L, namely in the continuum limit, the Wilson action becomes equivalent 
to the heat-kernel action ||I3| , p!^ . To be more precise, if we consider the action (P) 



with /?s = 0, in the large L and hence large f3t limit we have 

exp ( E E E ReTrt/^(f , t)V{x + /i, t)?7^(x, t + l)V\x, t) | 

\ S t=l/i=l / 

n (E^r-Xr(f/M(5^,t)n5^+A,t)f/J(5^,t+l)V^^(5',t))expf-;^')| ,(45) 



2N(3 



where r labels the irreducible representations of SU{N), their dimensions, 
their quadratic Casimir and Xr{U) the corresponding character of U . Consistently 
with previous notations we denote f3t by {3 in (|^) and following equations. The 



heat-kernel action at the r.h.s. of (^) was studied in detail in ref.s |0] and |l2 
which papers we refer the reader to for further details. The integration over the 
space-like links can be done exactly by using well known properties of characters, 
leading to the following action for the Polyakov loop: 



exp(5poi) = n E Xr{V{x + fi))xr{V\x)) exp f-^^ j 



(46) 
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where V{x) = Ylt=i^i^j'^) is the untraced Polyakov loop. By assuming as usual 
that the master field is invariant under translations, it was shown in [Q] that from 
one obtains the following partition function: 



N 



i=l 



exp 



/? , d-1 

—a 

L 12 



1 d 



(47) 



The integration variables 9i in eq. (^^ are the invariant angles of the Polyakov loop, 
and higher order terms in 9i in the exponent have been neglected. The partition 
function given in ( ^Tj) is again a Kazakov-Migdal model with a quadratic potential, 
and it admits the same semicircular distribution ( ^2|) as the partition function (|38|) 
obtained from our one plaquette model, but with (3h replaced by [3/L. This is in 
complete agreement with the leading order in L of the rescaling law given by eq. 
(p9|). The only approximation involved in deriving eq. ( ^7]) is the truncation of the 
power expansion in 9i at the exponent to the quadratic terms. In conclusion, the 
twisted one plaquette model and the heat-kernel model have in the large L limit the 
same solution for the eigenvalue distribution of the Polyakov loop, and the same 
critical value of the coupling constant. This confirms that the local extrema which 
were neglected in the saddle point evaluation of (^4|) are indeed irrelevant, at least 
in the large L limit and up to the critical value of the coupling constant. 

The relation between our one plaquette model and the heat-kernel model dis- 
cussed in IQ also provides a new and a priori independent way of estimating the 
critical coupling. In fact it was shown in ||^, and it can be recognised in (^6Dp|, that 
in the heat-kernel model the building block of the effective action for the Polyakov 
loop is given by the action of QCD2 on a cylinder of area inversely proportional to 
/3/L. So the critical value of Ph ^ P/L should correspond to the one where the 
Douglas-Kazakov phase transition [|T^ occurs. The point of the Douglas-Kazakov 
phase transition is known exactly as a function of the area of the cylinder and of 
the radius of the semicircular distributions at its boundaries. The radius of the 
distribution is obtained from the solution of (0) and it is given, as before, by (^31). 
The critical coupling is determined now by the radius where the Douglas-Kazakov 
phase transition occurs (see eqs. (33) and (42) in [Q]) rather than by looking at the 
point where the radius becomes complex. The results of Table I show that all the 
above methods give, at least for low d, values for the critical jSn which are very close 
to each other, and that the ones obtained from the analysis of the Douglas-Kazakov 
phase transition give consistently a better (higher) lower limit. 

It is well known that the Douglas-Kazakov phase transition is due to the con- 
densation of "instanton" contributions, namely to configurations where the initial 



^The expression under the product sign in ( pq ) is the partition function for QCD2 on on a 
cyhnder with boundary conditions given by the Polyakov loops in x and x + p.. 
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Table I: Critical value of Ph calculated in different ways: (a) limiting value for real 
radius, (b) limiting value for r < vr and (c) value corresponding to the DK phase 
transition in QCD2 §j. For d>3 the solution becomes unstable before the radius 
reaches the value vr. 







(a) 


(b) 


(c) 


d = 


2 


0.311 


0.314 


0.321 


d = 


3 


0.218 


0.218 


0.226 


d = 


4 


0.184 


* 


0.192 


d = 


5 


0.167 


* 


0.173 



and final value of one or more eigenvalues of V differ of more than 2tt. We argue 
that there is a precise correspondence between the instantons in QCD2 on a cylin- 
der and the classical solutions of the one plaquette model (extrema of the integral 
I{0,9 — ^)) which were neglected in the saddle point calculation of the previous 
section. These correspond to permutations which are not of the type of eq. (pBD, 
because one or more eigenvalues of V are not mapped from a bunch a to a + 1. The 
argument goes as follows: in the large L limit the eigenvalue distribution of our 
solution can be represented on an infinite line as a sequence of bunches, each made 
of n = N/L eigenvalues^, at intervals of 27r. Above the phase transition different 
bunches do not communicate with each other, except for the canonical shift of 1 
due to the twist, and the partition function is just the one corresponding to a single 
bunch raised to the power L (see eq. ([38| ) ). The picture is quite similar for the 
eigenvalues of the Polyakov loop in the heat-kernel model. Their distribution can 
be represented on an infinite line as a periodic distribution of period 2tt, that before 
the phase transition consists of separate bunches of eigenvalues. Let us compare 
now the contribution of classical solutions of the instanton type in the two models. 
Instanton contributions to the partition function of QCD2 on a cylinder can be 
calculated from the following equation (see eq. (8) in |T^) 



oc 



Y.Xr{-6)Xr{0) exp 

r 

^{l-N)/2 

— exp 



tCr 



1 



N 



i=l 



2nsif 



(48) 



where t can be identified with 



by comparing eq. (|48|) with (^Bf). The integers 
Si correspond to the winding numbers of the eigenvalues, so the contribution of an 
instanton with winding number s is proportional to exp(— 27r^/9A^s^/L). 



^We refer here to the eigenvalues normahsed as fa in (p 
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Let us compare this result with the non perturbative contribution (|29| ) corre- 
sponding to an eigenvalue being mapped by a permutation P from a bunch a to 
the bunch a + s + 1 in the twisted one plaquette model. The r.h.s. of (^) coincide 
with the instanton contribution calculated above, if one takes into account that the 
number of eigenvalues is here n = N/ L instead of N. 

We can conclude that the mechanism underlying the phase transition is the 
same in QCD2 on a cylinder and in the twisted one plaquette model. In particular 
in the latter model only fluctuations around the twist eating configurations need to 
be considered for /3 larger than its critical value; when the critical value is attained 
a condensation of local extrema of the action takes place, that is responsible for the 
phase transition studied in this section. 

4.2 Quantum rescaling 

Let us consider again the rescaling law given by eq. (0), and write it as: 

P{L) = Lf3H + \ , (49) 

where in the l.h.s. we replaced (3 with to remark the fact that the coupling 
refers to a model with L time-like links. On the other hand we have L = prit (see 
footnote in section 2.2), where p is the asymmetry parameter and the number of 
links in the time direction in the equivalent symmetric lattice. By writing eq. (|49|) 
with L = prit and with L = rit, and by eliminating jSn from the two equations, we 
find 

(3(pn,) = p^f3{nt)-\^ + ^ . (50) 

This equation gives the riscaling of the coupling constant induced by varying the 
asymmetry parameter p, and it should be compared with the first of eqs. (^ where 
/5symm ^ Pint) and Pt{p) — Pipnt)- The function Cr{p) of eq. (|^) was calculated 
by Karsch in |0. Its behaviour for large p is of the type 



ApCr{p)=a'ip + al + o(^^'^ 



(51) 



where the values of the coefficients and al in the large- iV limit can be calculated 
from the eq. 2.25 of 

a° = -0.2609 ; «t = J • (52) 



We can extract for comparison the same quantities from the rescaling ([50|), following 
from our effective model. We find : 



^t = -\ ; "t = i • (53) 
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The agreement is quite remarkable, considering that Karsch's calculation took into 
account the space-like plaquettes, and it confirms a feature that already emerged 
in the case of SU(2) 0], namely that the corrections to Cr{p) due to the space-like 
plaquettes are relatively small. The agreement between ( p3D and (|52D also represents 
a further check that only the fluctuations around the classical vacuum are relevant 
in the deconfined phase. 



4.3 Evidence for unbroken Zl symmetry 



In section 4 we solved the model described by the action (pq ) by assuming that the 
solution is translational invariant in the one dimensional space labelled by the index 
a in ipa,a- In other words we assumed that the symmetry of the vacuum (^) is 
not broken by the quantum fluctuations. Although this is justified in the large L 
limit, it would be desirable to have an independent proof valid for all values of L. 
In this section we shall prove that for any L the partition function (|36D becomes 
singular at the critical value of j3 derived in the previous section, namely at the value 
of /3 for which the radius of the large L semicircular solution becomes complex. 
To begin with, let us notice that the partition function ( p6|) can be obtained 



by integration over the unitary matrices U of the following model of Kazakov- 
Migdal type: 



a=l 



exp{-n/3^Tr [j] - J] $,f/W$,+i[/J'^) t] } 



a, 1^1 



(54) 

where the matrices are normalised so that their eigenvalues are Lipa,a, and the 
parameters Ph and are given by eqs. (|39|) and (Pl). 

The structure of the model (0) is illustrated in Fig. 2; it looks like a KM model 
defined on a compact one-dimensional lattice of length L, but with d links joining 
two neighbouring sites. 



If one integrates eq. (54) over the hermitean matrices $a rather than over the 
unitary matrices U, one obtains the "induced gauge model" [|TU| 



z./[nn«.r]exp{ii:J|^} 



1 ^ |Tr;7r 
/r(m,2)'r 



(55) 



where F is a closed path on the lattice of Fig. 2 and /r its length. Ur is the ordered 
product of the link matrices Ujf"^ on the path F. 

Consider now the quantity = ^(Tr$^), which is a measure of the width of 
the a'th bunch and is proportional, in case of a semicircular distribution, to the 
squared radius r^. It is easy to show, by using eq. (|5l) , that 



pI 



1 



n 



(Tr$^) 



■(E 



|Tr[/rP 



2Vr + 



T ) 



(56) 
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Fig. 2: Structure of the effective KM model corresponding to tlie case of L generic. 
Each hnk ^» — > — > V"^ represents a term ^aUjf''^ ^a+iUji^'^^ in the action. 



where the sum is over all paths F beginning and ending in a and the last expectation 
value is taken with respect to the partition function (^51). In the large-n limit we 
can use the factorisation property: (|Tr?7rP) = (Trf/r) (Tr?7p). Moreover due to 
the following symmetry of the action (^): 

f^i"^^ - Ci"^f^"\ Ci'^^ez„ , (57) 

which is unbroken in the weak coupling (large (3h) regime, we have: 

f if r is not a "tree" , , 

(Trf/r) = <^ (58) 
I Tr 1 = n if r is a "tree" , 

where by "tree" we denote any closed path T that contains Ujf'^ and Ujf^'^ the same 
number of times for any fi and a. By using this property, eq. (0) reduces simply 

1 (/) 

Pa ~ IT' X/ 7 ' '-^^^ 

Ph I (m^)'+i 

where n{l) is the number of trees of fixed length / on the lattice of Fig. 2 starting 
end ending in the point a. The instability of the weak coupling configuration occurs 
at the value of such that the series at the r.h.s. of eq. (|59D diverges. 

The basic point now is that^ irrespectively of the choice of the lattice on which 
the Kazakov-Migdal model (with quadratic potential) is defined, the expectation 
value of (Tr<l>^) on a site of the lattice is given by (^) where n{l) is the number of 
trees of fixed length / on the lattice we considered, starting end ending in the chosen 
point. On the other hand, n{l) depends only from the coordination number of the 

®For related discussions, see 1lS|l. 
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lattice, and hence it is the same for the lattice of Fig. 2 and for a (i-dimensional 
hypercubic lattice. In the latter case we know the exact solution, a semicircular 
distribution of radius given by (|^); (Tr$^) is proportional to the square of such 
radius and it becomes singular at the value of jSn given in eq. (^l]). In conclusion, 
although we have not been able to solve explicitly for all L's the model (|36| ) without 
making the explicit assumption that the symmetry of the vacuum is unbroken, 
we have proved that (^6]) becomes singular at exactly the same critical value of Ph 
as our Zl invariant solution. This shows that no phase transition associated to the 
breaking of the Zl symmetry occurs and it is a conclusive evidence that the Zl 
invariance of the vacuum is preserved. 

5 Conclusions 

We have considered the one plaquette matrix model corresponding to a hot twisted 
Eguchi-Kawai model in which the time-like plaquettes are neglected, and we have 
found a solution by calculating the quadratic quantum fluctuations around the 
classical vacuum, namely around the twist eating configurations. The solution for 
the time-like link variables consists of L bunches of eigenvalues, centered around 
the L^^ roots of unity, each bunch with a semicircular distribution of eigenvalues. 
This eigenvalue distribution, although rather unusual, does not come as a surprise 
in this context because it corresponds to a spreading of the eigenvalue distribution 
of the classical vacuum of the Twisted Eguchi-Kawai model. Our solution becomes 
unstable at a critical value of the coupling constant, where a phase transition occurs. 
The crucial point here is that our calculation neglects the contributions of the other 
extrema of the action, so that we do not know, without further information, the 
range of validity of the solution. This problem is solved by observing that in the 
hamiltonian limit {L oo) the original model is mapped into one with a heat-kernel 
action, which was solved in a previous paper 0]. As the two solution coincide we 
can safely say that the contribution from the relative extrema of the action can be 
neglected until the instability sets in. We can actually infer that the instability of 
the solution is determined by a condensation of the local extrema of the action. 
In fact it was shown in that the phase transition in the heat-kernel action 
is mathematically the same as the Douglas-Kazakov phase transition in QCD2, 
which is known to be due to a condensation of classical solutions (instantons). An 
interesting possibility is that the full theory, namely the one including the space-like 
plaquettes, might follow exactly the same pattern. It has already been suggested 
||20[| that a condensation of the extrema of the classical action is responsible for 
confinement in the Twisted Eguchi-Kawai model; the results of the present paper 
suggest that the fluctuations around the twist eating configurations are the only 
ones that need to be taken into account above the deconfinement transition and 
that they determine the master field of the theory in that phase. An instability of 
the corresponding solution would then give a lower bound for the critical coupling. 
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without having to take into account exphcitly the comphcated structure of the 
local extrema, whose condensation is ultimately responsible of the deconfinement 
transition at the critical temperature. The comparison between column (a) and (c) 
of Table I shows that the instability of the solution in the deconfined phase and the 
condensation of "instanton" solutions are indeed very close to each other, at least 
in our simplified model. 
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A Appendix 



In this Appendix we derive the summation formula given in eq. which is 

crucial in the computation of the quadratic effective model of section ||. 
This formula can be derived starting from the identity 



_[_[ 2sin(x 

1=0 



L 



2 sin Lx 



(60) 



which in turn can be proved by expanding the sine functions in exponentials and 
by noticing that the L^^ roots of unity satisfy the property: 

0, n = 1,2,...,L-1. . (61) 



k=0 



By taking the logarithm of both sides of eq. (|60D one obtains the relation: 



L-l 

1=1 



2 sin(a; + 



log 



sin Lx 



smx 



(62) 



which is the "generating identity" of a set of non-trivial trigonometric summation 
formulas. In fact by expanding the two sides of eq. (^) in powers of x we have: 



0(x) 



Oix" 



Ef=i' cot f 







(63) 



E 



1=1 



sm 



-2 111 



where the last equation coincides with (|33|). It is worth to remark here that, by 
taking appropriate combinations of higher order identities, one can evaluate from 
the power expansion of ( p^ ) all the sums of the form Y,i=i sin~^'^ ^. For instance, 
form the fourth and sixth order in x one gets: 



L-l 



^sm — 

1=1 ^ 



{L^ - 1) 2(L2 - 1) 



45 



(64) 



21 



and 

References 

T. Eguchi and H. Kawai, Phys. Rev. Lett. 48 (1982) 1063. 

For a review on Eguchi-Kawai models see S.R. Das, Rev. of Modern Physics 
59 (1987) 235 and references therein. 

M. Billo, M. Caselle, A. D'Adda, L. Magnea and S. Panzeri, Nucl. Phys. B435 
(1995) 172. 

M. Billo, M. Caselle, A. D'Adda, and S. Panzeri, Nucl. Phys. B472 (1996) 
163. 

F. Karsch, Nucl. Phys. B205 (1982) 285. 

F. Klinkhammer and P. van Baal, Nucl. Phys. B237 (1984) 274. 

K. Fabricius and C. P. Korthals-Altes, Nucl. Phys. B240 (1984) 237. 

A. Migdal, Sov. Phys. JETP 42 (1976) 413, 743. L.P. Kadanoff, Ann. Phys. 
(N.Y.) 100, (1976) 359, and Rev. Mod. Phys. 49 (1977) 267. 

I.I. Kogan et al, Nucl. Phys. B395 (1993) 547. 

V.A. Kazakov and A.A. Migdal, Nucl. Phys. B397 (1993) 214. 

D. Gross, Phys. Lett. B293 (1992) 181. 

K. Zarembo, Mod. Phys. Lett. AlO (1995) 677. 

J.M.Drouffe, Phys. Rev. D18 (1978) 1174. 

P. Menotti and E. Onofri, Nucl. Phys. B190 (1981) 288. 

M.R. Douglas and V.A. Kazakov, Phys. Lett. B319 (1993) 219. 



M. Caselle, A. D'Adda, L. Magnea and S. Panzeri, |iep-th/9309107| , Proc 



Trieste Summer School in high energy physics and cosmology. World Scientific; 
D. Gross and A. Matytsin, Nucl. Phys. B437 (1994) 541. 

Harish-Chandra, Am. I. Math. 79 (1957) 87; C. Itzykson and J.B. Zuber, J. 
Math. Phys. 21 (1980) 411; M.L. Mehta, Comm. Math. Phys 79 (1981) 327. 

M. Caselle, A. D'Adda and S. Panzeri, Phys. Lett. B293 (1992) 161. 



22 



[19] Y. Makeenko and H.W. Pe, ITEP-TH-13/95, |hep-th/960ll39| ; Y. Makeenko, 
Mod. Phys. Lett A8, No. 3 (1993) 209. 

[20] P. van Baal, Comm. Math. Phys. 9 (1983) 1. 



23 



